% Script calculates the solution of integral equations for the current
% distribution and input impedance of the antenna of radius a located on
% the surface of plasma column of radius b in the case of non-resonant
% plasma inside the column

clear all;
clc;
close all;
format short e;

% set constants
CONSTS = [];
% f
% (18e+6 Hz --- 400e+6 Hz @ B_0 = 6)
% (30e+6 Hz --- 1200e+6 Hz @ B_0 = 10)
CONSTS.f = 1050e+6; % 
CONSTS.c = 3e+10; % cm/s
CONSTS.a = 2.5; % cm
CONSTS.gamma = 0.577;
CONSTS.e = -4.80320427e-10;
CONSTS.m_e = 9.10660e-28;
CONSTS.B_0 = 200; % 620 200 5 --- 800 gauss 200
CONSTS.N = 1e+11; % 1e+11 --- 1e+13 cm^-3
CONSTS.d = (CONSTS.a)/50; % 0.1 cm --- 0,01 mm
CONSTS.delta = (CONSTS.a)/50; % 4.39e-2 rad for a = 0.4 --- 1.9 cm
CONSTS.phi_0 = 0;
CONSTS.eta_a = 1;
CONSTS.use_S = false;

CONSTS = update_CONSTS(CONSTS);

%--------------plot the current distribution and phase---------------------
printout_of_amplitude_and_phase(CONSTS, false, true);


